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Abstract. We make predictions for the kaon interferometry measurements in 
Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC). A first order 
phase transition from a thermalized Quark-Gluon-Plasma (QGP) to a gas of 
hadrons is assumed for the transport calculations. The fraction of kaons that 
are directly emitted from the phase boundary is considerably enhanced at large 
transverse momenta Kj* ~ lGeV/c. In this kinematic region, the sensitivity 
of the Ro\it/R a ide ratio to the QGP-properties is enlarged. The results of the 
1-dimensional correlation analysis are presented. The extracted interferometry 
radii, depending on ify, are not unusually large and are strongly affected by 
finite momentum resolution effects. 



1. Introduction 

We discuss predictions for the kaon interferometry measurements in Au+Au collisions 
at the Relativistic Heavy Ion Collider (RHIC) that accelerates the nuclei up to nucleon- 
nucleon center-of-mass energies of \/s NN = 200 GeV. Correlations of identical particle 
pairs, sometimes also called HBT interferometry, provide important information on 
the space-time extension of the particle emitting source as produced for example 
in ultrarelativistic heavy ion collisions [Q. In this case, QCD lattice calculations 
have predicted a transition from quark-gluon matter to hadronic matter at high 
temperatures. For a first-order phase transition, large hadronization times have been 
expected due to the associated large latent heat as compared to a purely hadronic 
scenario 0, ||, Q . Entropy has to be conserved while the number of degrees of freedom 
is reduced throughout the phase transition. Thus, one has expected a considerable 
jump in the magnitude of the interferometry radius parameters and the emission 
duration once the energy density is large enough to produce quark-gluon matter. Two 
alternative scenarios of the space-time evolution, with and without a phase transition, 
are illustrated in Fig. 1 in the z-i-diagram. After the collision of the two nuclei, each 
with nucleon number A, the system is formed at some eigen-time r (indicated by the 
hyperbola) and the initial expansion proceeds either in a hadronic state (left-hand side) 
or in a state dominated by partonic degrees of freedom, for example a quark-gluon 
plasma (QGP) (right-hand side). In the latter case, the formation of a mixed phase, 
leads to large hadronization times and thus to rather long emission durations. The 
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freeze-out is denned as the decoupling of the particles, i.e., the space-time coordinates 
of their last (strong) interactions. Summarizing this idea illustrated in Fig. 1, the 
interferometry of identical particle pairs and in particular the excitation function of 
the interferomtry parameters have been considered as an ideal tool to detect the 
existence and the properties of a transition from a thermalized quark-gluon plasma to 
hadrons. Quantities of great interest are for example the critical temperature T c or 
the latent heat of the phase transition. The prolonged emission duration should lead 
in particular to an increase of the effective source size in the outward direction, i.e., 
parallel to the transverse pair velocity. One also expects correlation lengths depending 
on the specific entropy of the collision. For recent reviews see for example J|, ||]. 



N, K, II 




Figure 1: Illustration of the space- 
time evolution in the z-i-diagram 
with (right) and without (left) 
a phase transition. Proceeding 
through a first-order phase transi- 
tion with a large latent heat should 
lead to large hadronization times, 
thus yielding eventually large in- 
terferometry radius parameters and 
emission duration. 



2. Space-time evolution and kaon interferometry 

Here, we discuss calculations based on a relativistic two-phase dynamical transport 
model that describes the early quark-gluon plasma phase by hydrodynamics and the 
later stages, after hadronization from the phase boundary of the mixed phase, by 
microscopic transport of the hadrons. For the initial dense (hydrodynamical) phase 
of a QGP a bag model equation of state exhibiting a first order phase transition 
is employed j3| [?]]. Hence, a phase transition in local equilibrium that proceeds 
through the formation of a mixed phase, is considered. A crossover [||, i\ or a rapid 
out-of-equilibrium phase transition similar to spinodal decomposition [9] may yield 
smaller radii and emission times. In the hadronic phase, resonance (de)excitations 
and binary collisions are modeled based on cross sections and resonance properties 
as measured in vacuum Jh| |ll). A detailed description of this relativistic hybrid 
transport model and its predictions can be found elsewhere |§, |lf , • Th e 

particular 



differences and advantages compared to pure hydrodynamical calculations are for 
example due to the explicit calculation of the freeze-out phase-space distributions. 
Moreover, the system evolution in the later dilute stages is based on cross sections 
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and resonance (de)excitations are possible. The assumption of ideal hydrodynamics 
might be not valid anymore in these later stages close to a freeze-out criterion like 
a hxed temperature Tf r . Another important example will be shown below (Fig. 5), 
demonstrating the possibility of direct emission from the early phase in this transport 
model used here. The mean freeze-out times are strongly depending on the transverse 
momentum of the particles as illustrated in Fig. 2 for kaons and anti-kaons. Obviously, 
large transverse momentum kaons are correlated to early mean freeze-out times. The 
difference between kaons and anti-kaons is almost negligible, only at very small pt anti- 
kaons may be emitted a tiny instant later due a small finite nct-baryon number in the 
system. (The K~N cross section is considerably larger than the K + N cross section.) 
Otherwise the Kir cross section, dominated by the K*, yields almost identical freeze- 
out distributions for kaons and anti-kaons. Strangeness distillation []l5| due to local 
large net-baryon numbers leading to emission asymetries that are observable in K + K~ 
correlations |L6) is not considered here. 
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Figure 2: Mean freeze-out time tf 
as a function of transverse momen- 
tum p t for kaons and anti-kaons in 
Au+Au collisions at RHIC. High 
Pt kaons are correlated with early 
freeze-out times. 
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The coordinate system for the correlation analysis is defined by the long axis (z) 
(parallel to the beam axis), the out direction (parallel to the transverse momentum 
vector Kt = (pti + PT2)/2 of the pair), and the side direction (perpendicular to 
both). Due to the definition of the out and side direction, i? ou t probes the spatial 
and temporal extension of the source while i? s ide only probes the spatial extension. It 
has been suggested that the ratio R ut /-Rsido should increase strongly once the initial 
entropy density Si becomes substantially larger than that of the hadronic gas at T c Q . 
The Gaussian radius parameters can be obtained from a saddle-point integration over 
the classical phase space distribution of the hadrons at freeze-out (points of their 
last (strong) interaction) that is identified with the Wigner density of the source, 

s(a)@|@|. 

> 2 (K T ) Hy 2 )(K T ), 



-Rsidel^Tj 
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being the space-time coordinates relative to the momentum dependent effective source 
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centers. The average in (§)-(§) is taken over the emission function, i.e., 
, w % fd 4 xf(x)S(x,K) 

with K = {E K , K). In the osl system = (/3 t ,0,/3z), where (3 — K/E K and E K = 
\fm 2 + K 2 (on-shell approximation). For small i-t correlations, i.e. in particular at 
small Kt, Rout is increased relative to i? s ide if the duration of emission At = ^{t 2 ) 
is large || . Strong (positive) i-i-correlations or large spatial anisotropies in out- 
and side-direction {(y 2 ) > (J 2 )) may, in principle, lead to i? ut < -Rsidc- This is not 
seen within the model scenario discussed here. 

From eqs. (1) and (2) one can easily calculate the ratio R u_t/R s idc for kaons. 
Kaons, compared to pions, are expected to be less contaminated by resonance 
decays jl8|, |oL In addition the kaon phase space density is much smaller than the 
pion density | 2l| . As a consequence, higher multiparticle correlation effects (that 
may be important for pions at RHIC energies ^2|) are possibly under much better 
control for kaons. Fig. 3 shows the kaon i? ut/-Rsido ratio for SPS (initial conditions 
Ti = 1 fm/c (thermalization eigen-time) and s/ps =45 (specific entropy density)) and 
RHIC (n = 0.6, s/p B = 200) @|. The bag parameter B is varied from 380MeV/fm 3 
to 720MeV/fm 3 . This corresponds to a variation of the critical temperature from 



160 MeV to T, 
2.4 



200 MeV. 



2.2 
2.0 

a 1.8 



• Ro/R s T=160 RHIC 

□ Ro/R s T=160 SPS 

A RJR S T=200 RHIC 

Ro/R s T=200 SPS 




Figure 3: Rout/R B iac as obtained 
from eqs. (1) and (2) for kaons at 
RHIC (Au+Au, yfi NN = 200 GeV) 
(full symbols) and at SPS (Pb+Pb, 
y/s NN — 17. 4 GeV) (open sym- 
bols), as a function of Kt for criti- 
cal temperatures T c ~ 160 MeV and 
T c ~ 200 MeV, respectively. 
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For small Kx, the ratio R QU t/Rside increases slowly compared to the faster rise 
for pions p3] j. This is due to the different masses yielding different flow velocities 
at the same Kt (see eq. (2)). Apparently, the sensitivity to the critical temperature 
T c and the specific entropy density increases strongly with Kt- For higher T c the 
hadronization is faster but the subsequent hadronic rescattering phase lasts longer. 
It was shown that this dissipative hadronic phase dominates the radii for pions p3| . 
The ratio R ut/Rside is always larger than unity and reaches values on the order of 
« 1.5 — 2 at large Kt- This value is considerably smaller than former expectations (e.g. 
[0. ^fj). On the other hand, first RHIC data for pion correlations show ratios that do 
not increase with Kt and are even smaller than unity ]25|, p6| . This completely new 
behaviour has not been seen at SPS energies (see e.g. p7|). However, new data by 
the CERES collaboration show a similar trend (28) . This observation would hint at a 
rather explosive scenario with very short emission times, not compatible with a picture 
of a thermalized quark-gluon plasma hadronizing via a first-order phase transition to 
an interacting hadron gas |23||. Rather a shell- like emission as illustrated in Fig. 4 
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would be prefered. Thus, the further study of two-particle interferometry will provide 
extremly important information e.g. on the hadronization process or the question of 
thermalization in ultrarelativistic heavy ion collisions. Also, the comparison of pion 
and kaon interferometry data may further clarify the situation. 




Observer 



Figure 4: Illustration of a shell- 
like emission. The surface emission 
geometry corresponds to small values 
of the ratio Rout/Rsidc, indicated 
by the two dashed lines in the 
emission volume element relevant for 
an observer. The dashed line in the 
direction to the observer corresponds 
to the out direction, the orthogonal 
line is the side direction. 



"Source" 

Kaons are less distorted by long-lived resonances and escape the opaque hadronic 
phase easier (compared to pions). About 30% of the kaons at Kt ~ 1 GeV/c are 
directly emitted from the phase-boundary (for T c w 160 MeV) as can be seen in Fig. 5. 
Complementary, we have seen already that large Kt kaons and their i? ou t / i? s ide ratio 
exhibit a strong sensitivity on the QCD equation of state. The fraction of resonance 
decays (K*'s,. . .) is still quite large and decreases with Kt from about 70% to 50%. 
For the higher T c (~ 200 MeV) , hadronization is earlier and the hadronic phase lasts 
longer, such that the system is more opaque for direct emission than in the lower T c 
case. 
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Figure 5: Fraction of kaons r rigm that origin from a particular reaction channel 
prior to freeze-out. These are resonance decays (full circles), direct emission from the 
phase boundary (stars), elastic meson-meson (diamonds), or elastic meson-baryon 
(open circles) collisions. The upper and lower diagrams are for RHIC and SPS initial 
conditions for T c ~ 160 MeV (left) and T c ~ 200 MeV (right), respectively. 
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We now calculate the one-dimensional kaon kaon correlation functions (^(Qmv) 
for various if^-bins and determine the corresponding fit parameters. The 
corresponding results for the three-dimensional analysis in the out-side-long coordinate 
system are presented in [ fTi] , The parameters Ri nv and Ai of the correlation functions 
are obtained by fitting a Gaussian as 

C 2 {Q im ) = 1 + Ai exp(-g2 nv i? i 2 nv ) . (6) 

The correlation functions themselves are calculated from the phase space distributions 
of kaons at freeze-out using the correlation after burner by Pratt g [TJ). The 
correlation functions can be calculated as 

/ d 4 xS(x, K) / d 4 yS(y, K) exp (2ik ■ (x - y)) 



C 2 (p 1)P2 )^l+ |/ d ^0r,K)|2 

with 2K = pi + P2, 2k = pi — p 2 , and 2k° = E Pl — E P2 . 

10 i — ■ — . — ■ — . — ■ — . — 



(7) 
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Figure 6: Parameters _Ri nv and Ai (as obtained from a fit according to eq. (6)) as a 
function of Kt for kaons in Au+Au collisions at RHIC for critical temperatures T c ~ 
160 MeV (circles, squares) and T c ~ 200 MeV (triangles, diamonds), respectively. 
Calculations with and without taking momentum resolution effects into account are 
shown. 

The results for RHIC are shown in Fig. 6. The Ai intercept parameters 
are multiplied by a factor 10. The correlation functions have been calculated for 
T c w 160 MeV (circles and squares) and T c w 200 MeV (triangles and diamonds). 
Moreover, calculations with and without taking finite momentum resolution (f.m.r.) 
effects have been performed. In the first case, the true particle momentum p obtains an 
additional random component. This random component is assumed to be Gaussian 
with a width dp. A f.m.r. of w 2% of the center of each Kt bin is considered, a 
value assumed for the STAR detector [^5] . The relative momenta of pairs are then 
calculated from these modified momenta. However, the correlator is calculated with 
the true relative momentum. i?i nv is on the order of 5.5 — 6.5 fm without and 4.5 — 6 fin 
with taking f.m.r. into account. Hence, -Rj nv is reduced by the f.m.r.. This reduction 
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grows with Kt- For the one-dimensional fit parameters, there is hardly a sensitivity 
to the critical temperature. However, the known trend (larger radii for higher T c due 
to a prolonged hadronic phase) is also indicated in these results. The sensitivity to 
the specific entropy density, another initial condition, is known to be rather weak | g3| |, 
even for a change from SPS to RHIC initial conditions, due to the dominance of the 
late soft hadronic interactions. Hence, the differences in the correlation parameters 
between, for example, the RHIC 130 GeV and 200 GeV energies are minute; (the 
number of charged particles at midrapidity increases only by about 15%). Preliminary 
analysis of STAR KK data (not yet corrected for f.m.r.) at \/s NN = 130 GeV gives a 
value of R inv = 4.5 ± 0.3 fm (at 150MeV< K T < 400 MeV and rather central (11%) 
collisions) (Ai = 0.92 ±0.13) 1^91 , which is almost compatible to the calculations that 
take f.m.r. into account. The Ai intercept parameters are almost constant without 
f.m.r. effects but strongly decrease when taking them into account. For illustration, 
let us consider particle pairs with a small true relative momentum q that will be 
redistributed by the random momentum smearing to, on the average, larger relative 
momenta q. Thus, the area around q ~ is depleted from these true low momentum 
pairs that carry the full correlation strength. Complementary, some pairs with a 
larger true q (and thus with a weaker correlation strength) are redistributed to the 
small q area. As a consequence, the correlation strength is transported to larger q 
values. Thus, the correlation function gets broader and the radius parameters become 
smaller. The strong decrease of Ai with Kt is due to the absolute larger momentum 
smearing at high Kt- True correlated pairs at low q are transported more likely 
to larger q values. This reduced correlation strength in the first q bins causes the 
small Ai values at larger Kt- Recent experimental data (corrected for f.m.r.l) for 
central Pb+Pb collisions at </s = 17.4 GeV yield R inv = 6.2 ± 0.4 ± 0.9 fm at low 
transverse momenta ((p t ) = 250 MeV) and i?i nv = 3.3 ± 0.6 fm at high transverse 
momenta ((pt) = 910 MeV) ]3C[ ] . The corresponding Ai parameters are reported as 
Ai = 0.84 ±0.11 ±0.19 and Ai = 0.47 ±0.12 for the low and high p t set, respectively. 

3. Summary 

The calculation of kaon correlation parameters for Au±Au collisions at RHIC energies, 
assuming a first-order phase transition from a thermalized QGP to a gas of hadrons 
shows 

• an increasing i? ou t/-Rsido ratio with Kt (larger than unity), 

• an increasing sensitivity of the R OVL t/Rsidc ratio with Kt, 

• a strong (30%) direct emission component from the phase boundary at large 
K T ~ 1 GeV/c, 

• no unusually large radius parameters R- mv that are reduced due to finite 
momentum resolution effects, 

• Ai intercept parameters that get strongly reduced due to finite momentum 
resolution, 

• the reduction of i?j nv and Ai grows with Kt- 

The kaon interferometry measurements at RHIC (at high Kt) will be, in combination 
with the pion data, an utmost important and valuable probe of the space-time 
dynamics (close to the phase boundary) and to the properties of the phase transition. 
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